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Abstract: We calculate the jet shape and the jet cross section in heavy ion collisions using 
soft-collinear effective theory (SCET) and its extension with Glauber gluon interactions 
in the medium (SCETg). We use the previously developed framework to systematically 
resum the jet shape at next-to-leading logarithmic accuracy, and we consistently include 
the medium modification by incorporating the leading order medium-induced splitting 
functions. The calculation provides, for the first time, a quantitative understanding of the 
jet shape modification measurement in lead-lead collisions at = 2.76 TeV at the 

LHC. The inclusive jet suppression is also calculated within the same framework beyond 
the traditional concept of parton energy loss, and the dependence on the centrality, the 
jet radius and the jet kinematics is examined. In the end we present predictions for the 
anticipated jet shape and cross section measurements in lead-lead collisions at .^snn ~ 5.1 
TeV at the LHC. 






1. Introduction 


The Run II of the Large Hadron Collider (LHC) has started this year, and the nucleon- 
nucleon center of mass energy of the lead-lead collisions is planned to be boosted above 
5 TeV. More energetic jets will be produced, which make the LHC an excellent arena 
for the study of jet physics in heavy ion collisions. On the other hand, one of the top 
priorities of the heavy ion program at the LHC is to study the properties of the hot, dense 
medium that is produced during the collisions and referred to as the quark-gluon plasma 
(QGP). Many interesting questions are under investigation, for example, how the medium 
is created and subsequently thermalizes, what the dynamics of the medium is and how the 
medium properties evolve with energy and density. To answer these questions, one of the 
key insights is to study how the presence of the medium affects various hard processes as 
probes of the QGP. The medium responses to the perturbation from the hard probes may 
also contain useful information. 

With the abundance of jets and the high-resolution detectors at the LHC, one promis¬ 
ing approach to extract the medium properties is to study the modification of jets as they 
propagate through the medium. This is the jet quenching [T] phenomenum which has 
been observed at the Relativistic Heavy Ion Collider (RHIC) [2H7] and the LHC [8lt27j. 
Traditionally, jet quenching refers to the suppression of production cross sections for var¬ 
ious hadrons. The suppression was understood theoretically with several different ap¬ 
proaches [28H41| (for reviews with different perspectives, see [42H45j 1. They share the com¬ 
mon physical picture of the parton energy loss, but with emphasis on physics at different 
energy scales. It has been clear that hadron cross section suppression alone is not enough 
to identify and distinguish between different physical mechanisms in jet quenching |46II48| . 
and more differential and correlated measurements are needed [4911501 . 

The studies of jet substructure observables and production cross section can provide 
precision tests of the QCD dynamics. Historically, the studies of event shapes and cross 
sections in e'’'e“ collisions helped confirm the gauge theory structure of QCD |51H60| . The 
discovery of the asymptotic freedom of QCD at high energy allowed the understanding of 
jets using perturbative tools. Since then, accurate QCD calculations provide robust theory 
control of jet physics |61H63| . They have also allowed one of the most precise extraction 
of the strong coupling constant |64H68| . In heavy ion collisions, the intricate jet formation 
mechanism provides the sensitivity to the medium dynamics. It is expected that precision 
QCD calculations will be crucial in the understanding of jets, therefore allowing us to 
extract the medium properties reliably. 

In general, jet substructure observables are more sensitive to the details of the final 
state, jet-medium interactions. They allow us to disentangle the initial state, cold nuclear 
matter (CNM) effects and therefore are cleaner probes of the medium. On the other hand, 
different jet substructure observables are sensitive to radiation at different energy scales. 
By measuring, for example, from jet cross sections, jet shapes and jet fragmentation func¬ 
tions to angularities, jet masses and particle multiplicities, the in-medium jet formation 
mechanism across a wide range of energy scales can be examined. Jet substructure observ¬ 
ables and their medium modifications are also highly dependent on the partonic origin of 
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jets. Quark-initiated and gluon-initiated jets typically have different substructures because 
of the different color charges isgiio]. With different medium interactions, quark-jets and 
gluon-jets can even be considered independent probes. 

It should be noted that, jets are conventionally defined by a jet algorithm used in the jet 
reconstruction. An angular scale R is introduced, and jets are identified as exclusive objects 
in an event. While jet substructure observables probes directly the radiation inside jets, the 
jet cross section can give complimentary information about the radiation in the rest of the 
event. In heavy ion collisions, because of the huge underlying event background, a small 
radius is usually chosen in jet reconstruction. This leads to the significant suppression of the 
jet cross section because the radiation can easily be lost outside jets. Recent progress has 
been made on the precision calculations of jet cross sections for jets with small radii 
and a better understanding of the radius dependence of the cross section. 

Because of the sensitivity to physics at multiple energy scales, the theoretical calcu¬ 
lations of jet substructure observables are challenging tasks, which has been realized for 
a long time. At high energy, perturbative calculations are reliable, and the contributions 
from non-perturbative regimes are power-suppressed. However, because of the emergence 
of hierarchical scales, the calculations of jet substructure observables always suffer from 
large logarithms of the ratio between these scales. The resummation of the perturbative 
series thus becomes necessary, and the medium modification of jet substructure observables 
will have to be consistently included. 

We will study the jet shape and the jet cross section in heavy ion collisions. The jet 
shape m is one of the classic jet substructure observables and it probes the transverse 
energy profile inside a jet. Given a jet of size R reconstructed using a jet algorithm, the 
integral jet shape 'hj(r) is defined as the fraction of the transverse energy Et of the jet 
within a subcone of size r around the jet axis h, 

Tj(r) = . (1.1) 

2-^i, dif,<R 


Here din = \/{Vi ~ "HjetY + “ 4>jet)'^ is the Euclidean distance between the Tth particle 

in the jet and the jet axis on the pseudorapidity-azimuthal angle (ry, cp) plane. The average 
integral jet shape 'l'(r) = j{r) and its derivative, the differential jet shape 

p{r) = ^T(r), describe how the transverse energy inside jets is distributed in r. In heavy 
ion collisions, the modification of jet shapes is conventionally evaluated by the following 
ratio 


Mp(r) 


p^^{r) 
pPP{r) ■ 


( 1 . 2 ) 


For the modification of jet cross sections, we measure the nuclear modification factor Raa 
defined as follows. 


Raa 


d(j aa 

dpdpT 


(iVbi, 


da. 


pp 


dpdpT 


(1.3) 


where (A^bin) is the number of binary nucleon-nucleon collisions in an A-|-A reaction. 
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Historically, the jet shape was first resummed using the modified leading logarithmic 
approximation in m, and later in |75II76| using a different approach. In heavy ion collisions, 
the work [MIET] built upon [71] and studied the medium modification of jet shapes using 
the Gyulassy-Levai-Vitev formalism [341135] in the soft gluon limit. Other approaches 
include m, and Monte Carlo studies of jet shapes have also been performed [79ll8T| . On 
the other hand, calculations based on the framework of traditional perturbative parton 
energy loss [821485] have been able to describe the suppression of light-flavor jets and its 
radius dependence, heavy flavor jets [T9l[86l - [88] . as well as the enhanced asymmetry in di¬ 
jet and 7 /Z°-|-jet events observed in heavy ion collisions [T2l[2Tl[22l[2l] . Other approaches 
incorporate energy redistribution in the parton shower through various transport models, 
typically implemented in Monte Carlo simulations and event generators |89H96| . 

In |971[98], we initiated the studies of precision jet substructure and cross section 
calculations in heavy ion collisions using soft-collinear effective theory (SCET) |99H104| and 
its extension to include Glauber gluon interactions in the medium (SCEXq) [10511106] . The 
goal of such development is to incorporate recent advances in pQCD and SCET in accurate 
theoretical predictions of jet observables in heavy ion collisions. In this formalism, Glauber 
gluons are treated as background fields generated from the color charges in the QGP. They 
mediate the interactions between the collinear jet and the QCD medium. SCETq allows 
us to calculate the medium-induced splitting functions [10711108] beyond the traditional 
soft-gluon emission limit, which have been applied to describe the suppression of charged- 
hadron and neutral-pion cross sections at RHIC and the LHC [10911110] . In this paper, we 
will present the first perturbative calculation of the jet shape modification and jet cross 
section suppression beyond the traditional energy loss approach. 

The rest of the paper is organized as follows. In Sec. [7] we present the theoretical 
framework for the calculation of the medium-modified jet shapes and cross sections. We 
demonstrate how the medium-induced splitting functions are incorporated naturally in the 
jet energy function and the jet energy loss calculations, with the implementation of cold 
nuclear matter effects. In Sec. [3] we compare the theoretical calculations to the experimen¬ 
tal measurements in Pb-|-Pb collisions at = 2.76 TeV at the LHC. We also show 

predictions for the upcoming .^/snn ~ 5.1 TeV measurements. In Sec. 0] we summarize and 
give an outlook of future studies. 


2. Framework 

With the relevance of multiple energy scales in heavy ion jet physics, we aim to develop 
an effective field theory framework in which both the resummation of jet substructure 
observables, as well as the calculation of their medium modifications, can be performed 
consistently, and the precision can be systematically improvable. In this section we review 
the recently developed framework in SCET [9911103] for the resummation of jet shapes in 
proton collisions m- For the jet shape in heavy ion collisions, this framework is readily 
generalized [S5| to include medium modifications via Glauber gluon interactions [TTHirnH] . 
The medium-induced splitting functions [10711108] were calculated in SCEXq in a medium 
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model consisting of thermal quasi-particles undergoing longitudinal Bj or ken-expansions Q|. 
They aim at describing how the parton splitting in the full kinematic regime changes in a 
medium with the characteristic scale at the temperature T. The effect from a more subtle 
medium spacetime evolution is expected to be suppressed by the jet energy. 

SCET is an effective field theory of QCD suitable for describing processes producing 
jets with energetic collinear particles. In such processes, a dynamical scale mj emerges 
where mj is the jet mass and is much lower than the jet energy Q. SCET expands 
the full QCD contributions in the calculation with a systematic power counting. The 
power counting parameter A ~ TnjjQ is typically small and the leading-power contribution 
calculated in SCET is a good approximation. The power-suppressed corrections can also be 
included order by order systematically. In the calculation of certain types of jet substructure 
observables, an even lower scale rrijfQ emerges which is the scale of the ubiquitous soft 
radiation in an event. SCET separates physics at these hierarchical scales and makes 
the factorization of the physical cross section manifest. This allows the calculation to 
be factorized into pieces, each can be performed at its natural scale more easily. The 
contribution in each sector at any other scales can then be obtained by the renormalization 
group evolution. 

In SCET we have the following hierarchy of energy scales, 

Q ^ QX ^ QX^ > Aqcd- (2.1) 

In the presence of the medium at the characteristic scale T, the soft scale QX^ can be 
comparable to the medium temperature. For a jet along the direction h, in light-cone 
coordinates the momentum is expressed as p = {n ■ p,n ■ p,p±) where n = (l,n) and 
n = (1, —h). The momentum scalings of the collinear and soft particles in the jet, as well 
as the Glauber modes mediating the medium interaction, are 

p, = Q(l, A2, X),Ps = Q(X, a, X), pg = Q{X\ X\ A) . (2.2) 

Note that the Glauber modes are off-shell degrees of freedom describing the momentum 
transfer transverse to the collinear direction. They are not hnal state particles in the jet and 
therefore do not directly contribute to the measurement of jet substructure observables. 

^The medium-induced splitting functions depend on the density and the geometry of the medium. At 
the LHC the gluon density dominates and it is related to the soft, charged hadron rapidity density which 
is proportional to the number of participants in the optical Glauber model. The proportionality coefficient 
is extracted from the experimental measurements in lead-lead collisions up to ^snn = 2.76 TeV and 
extrapolated to 5.1 TeV. In particular, the gluon rapidity densities dN/dy we use in the calculation are 
as follows: at 2.76 TeV, we use 2640 (750) for central (mid-peripheral) collisions, while at 5.1 TeV we 
use 2750 (800). The medium-induced splitting functions are averaged within each centrality class over the 
spacetime position of the hard scattering process in the Glauber model, allowing all possible path lengths 
in the average. The event-by-event fluctuation of jet modifications and the correlation with the medium 
geometry will be studied in the future. 

^While many bulk properties of the medium require a complete hydrodynamic description, in this work 
we focus on jet physics above medium scales and evaluate the medium-induced splitting functions using the 
classic Bjorken-expansion model of the medium. Medium effects at lower energies on jets will be investigated 
in the future. 
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In the splitting of a collinear parton with momentum p = (a;, 0,0) into partons with 
momenta k = (xw, k‘j^/xoj, k^) and p — k, the formation time of such splitting process is 


x(l — x)uj 1 
tsplit ~ ~ 


(2.3) 


This is to be compared with the typical mean free path of the parton multiple scattering 
in the medium, 

tfree ~ ( 2 - 4 ) 

and tsplit can be comparable to tfree or even larger. This leads to the coherent multi¬ 
ple scattering of partons which underlies the Landau-Pomeranchuk-Migdal (LPM) effect. 
The induced radiation in the collinear direction is suppressed, resulting in the wide-angle 
bremsstrahlung m which is captured by the SCETq medium-induced splitting functions. 

In the context of the jet shape calculation, the integral jet shape has dominant contri¬ 
butions from the collinear sector in SCET, 


T(r) 


Er 


E^r + E^r 

E’k + E^ 


|l + 0(A), 

R 


(2.5) 


because E'^ = 0{Q) and E^ = 0{Q\). E^ and E^ are the contributions to the measured 
energies from collinear and soft radiation. The collinear nature of the jet energy leads to 
the simple factorized expression of the jet shape 193, 




Juj,Er{k) 
Juj,ER{k-) ’ 


( 2 . 6 ) 


which is the ratio between two jet energy functions, and p is the renormalization scale 
appearing in the calculation using dimensional regularization in d dimensions and the MS 
renormalization scheme. Here 


Ju,eM = Y,muo)\Xe){Xe\xum)Mxc) , ( 2 . 7 ) 

is the jet energy function which describes the average amount of jet energy inside the 
subcone of size r radiated from the parton with the collinear momentum h ■ p = oj. XtJ is 
the collinear jet field operator and Xc is the collinear sector constrained within the angular 
scale R used in the jet reconstruction. The measurement operator Er returns the energy 
inside the subcone of size r. More detailed information about the above expression can be 
found in 123, nevertheless it simply calculates the jet energy weighted by the probability 
using SCET matrix elements. 

In fact, the collinear sector in SCET is a boosted copy of QCD |104llll2HI14j . Realizing 
this, various next-to-next-to-leading order jet and beam functions were calculated bypassing 
the use of the more complicated collinear SCET Feynman rules [nsHns]. It was also 
pointed out that jet and beam functions can be calculated directly from the QCD splitting 
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Figure 1: Illustration of the phase space regions for the calculations of the jet energy function 
(left panel) and the jet energy lost (right panel) at leading order. In both cases, a collinear parton 
splits into partons with momenta k = {xPq , k^/xpQ , k±) and p — k. Depending on the kinematics 
of the splitting, the partons which can contribute to the energy measured inside the subcone of size 
r may change. Similarly, the partons reconstructed in the jet may also differ. 


functions |120j with proper phase space (PS) integrations. At leading order □, the jet 
energy function associated with parton i with the collinear momentum p splitting into 
k = {xuj, k‘j_/xu], A:_i_) and p — k can then be written as follows, 

JL,Erik-) = [ dx^^Vi^jk{x,k±)Er{x,k±) , ( 2 . 8 ) 

j,k J 

where Pi->.jkix, k±) are the collinear parton splitting functions and the d-dependence is 
suppressed. In the presence of the medium, 

kx) = v:%k{x, k±) + kx), (2.9) 

which is the sum of the vacuum Altarelli-Parisi splitting functions and their medium- 
induced modifications 0. Er{x,k±) is the measurement function associated with the jet 
energy function, 

E^[x, k±) = Adi -|- Ad2 + Ad3 -|- Ad4 , (2.10) 

where 

Adi = tan - — k±^ 0 ^a;(l — x) tan - — k^^ 0k,j. x p^ , (2-11) 

®While splitting functions and jet energy functions are calculated at fixed order, all-order contributions 
from vacuum parton splittings and multiple gluon emissions are included in the resummation at next-to- 
leading logarithmic accuracy. However, contributions from correlated multiple scatterings with the medium, 
i.e. higher-order terms in the opacity expansion, are not resummed in this work. They may affect the 
medium-induced splitting functions, and the impact on the jet shape modification will be investigated in 
the future. 

“^The medium-induced splitting functions in this paper were calculated using SCETg at the first order 
in the opacity expansion. 
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A ^2 = 0 {^x tan - — k±^ 0 (jt± — uj{l — x) tan 0^^ x , (2-12) 

A ^3 = 0 (^k± — wx tan 0^a;(l — x) tan - — k±^ 0^^ x (p^ — k^) , (2-13) 

Mi = 0^A;_l — cjx tan 0 ^A:_l — oj{l — x) tan 0^^ x 0 , (2-14) 

are the cases in which each particle is either inside the subcone of size r or not. 0krp is 

the phase space constraint corresponding to the reconstruction of jets of size R using the 
anti-lcT algorithm, 

0kT = 0— 3;) tan — — /c_l j . (2-15) 

The different phase space regions in the jet energy function calculation are illustrated in 
the left panel in Fig{T] It is then straightforward to calculate the medium modification of 
the jet energy function, 

(2-16) 

Note that, due to the LPM effect the medium modification in this treatment contributes 
as a power correction because V^ji,{x, A:_i_) is finite as k± —)• 0. The renormalization group 
evolution of the jet energy function is then the same as in vacuum, 


din fi 


C'iFcusp(Q^ s) In 


tan^ Y 


- 27* («^) 


^ i 


lo.E') 


(m)> 


( 2 . 17 ) 


where i = q,g with Cg = Cp and Cg = Ca the Casimir operators of the fundamental and 
adjoint representations in QCD. The anomalous dimensions Tcusp and 7* are given in m, 
and the solutions to the renormalization group evolutions have the same form, 

/^2 ^^^2 R \ 

= JL,Erif^jr)eW[-‘^CiS{pj^,p) +2Ai{pg^,p)] i -, ( 2 . 18 ) 


which describe how jet energy functions evolve from the jet energy scale pj^ to arbitrary 
renormalization scale p. Here S{i',p), Ai{v,p) and Ax^{v,p) are the renormalization-group 
evolution kernels in SCET and the explicit expressions can be found in m- The integral 
jet shape becomes 


K(r) = pjj^) + 2Ai{pj^, Pj^)] ( 

UJ,ER\k-jR) tan ^y 

As discussed in m. the same choice of jet energy scales 

T R 

Pj^ = oj tan - « Ej X r , pj^^ = uj tan — « Ej x R , 


( 2 . 19 ) 

( 2 . 20 ) 


where Ej is the jet energy, eliminates the large logarithms in the fixed order calculation 
of Er^k-jr) and E En^k^jn) because the medium modification does not introduce extra 
logarithms. The renormalization group evolution between pj^ and pjj^ then resum the 
logarithms Inpj^/pj^^ « Inr/R caused by the vacuum parton splitting. 
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With the resummation and the medium modification performed, the integral jet shape 
is the average over the jet production cross section. 


^'(r) 


*^total 


E 


IPS 


drjdpT 


drjdpT ^ 


(r) , 


( 2 . 21 ) 


with the proper phase space cuts (PS) on the jet px and p. With the necessity, we now 
turn to the calculation of the jet cross section in the medium. 

The jet cross section in nucleus collisions is calculated through the following expression 
with the Glauber modeling 


1 




(A^bin) dpdpT 


= E 

ijX ■ 


dxidX2ff'(xi; PCNM)ff(x2] fJ-CNM) 


do'ij^kx 
dxidx2dpdpT' 


( 2 . 22 ) 


where ff^{x; pcNu) are the effective parton distribution functions (PDFs) in the nucleus 
A, and pcNM as an argument in parton distribution functions parameterizes the cold nu¬ 
clear matter effect we implement, not to be confused with the factorization scale 0. The 
induced radiation of the initial-state and final-state partons can both cause a modification 
of the jet cross section. The former contributes as one of the cold nuclear matter (CNM) 
effects because the hot, dense medium is produced after the collision. It can give signif¬ 
icant contributions [123| to the attenuation of the jet production cross section in central 
d-l-Au collisions at RHIC |124| and p-|-Pb collisions at the LHC |125| . Although here we 
mainly discuss the production cross section of jets, measurements of the inclusive particle 
production in p-|-Pb collisions [126| and the quenching of hadrons with very high transverse 
momentum in Pb-|-Pb collisions at the LHC |110j are consistent with small CNM effects. 

The calculation of the CNM effect in this paper follows |122j in which the effective 
field theory techniques have not been applied yet. We implement the initial-state CNM 
effect as follows, 

fq{x', fJ-CNu) = j def^ fJ-CNu) , (2.23) 

/^cnm) = j def^ /^cnm) , (2.24) 

where ff{x) are the parton distribution functions of the proton i Pq,g(^', are the 

probability density for the initial-state quarks or gluons to lose a fraction e of their energies 
due to multiple soft gluon emissions as they propagate through a large nucleus. Note that 

®In this work we do not incorporate the intrinsic modifications of the parton distributions which are 
conventionally referred to as nuclear parton distribution functions (nPDFs). We only examine and refer to 
the cold nuclear matter effect from initial state parton energy loss mi- This is caused by the radiation 
induced in the presence of the nucleons in the nucleus by the RMS momentum exchange at the scale (J-cnm 
with the colliding partons. 

®Again we do not incorporate the nuclear parton distribution functions (nPDFs). As one example we 
focus on effects from initial state parton energy loss on jet shapes and cross sections. Distinguishing and 
combining various initial state effects in theoretical calculations, as well as detailed studies of |125ill27] 
where strong constraints on initial state effects can be made, are left for future work. 
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in this case the modifications of the effective PDF of the nucleus moving in the positive 
lightcone direction are caused by the interaction with the nucleus moving in the negative 
lightcone direction. The scale fJ-cNM, which will later show up in the Results section, is 
the scale of the RMS momentum exchange between the colliding partons and the nucleus 
(not to be mixed with the factorization scale in the PDFs). The parameter controls the 
strength of the CNM effect implemented in the calculations. 

The CNM effect can essentially be organized as the modification factors associated 
with the parton distribution functions, 

fJ-CNu) = fi{x) X i?j(x;/rcNM) ■ (2.25) 

Note that, because the PDFs typically die off as x —)• 1, the factors Ri are smaller than 
1 in this region due to the energy loss of initial state partons. This can cause partly the 
attenuation of the jet cross section in nucleus collisions. 

Similarly, because the jet cross section dies off as the jet transverse momentum pT 
increases, the jet energy loss due to the final-state, medium-induced radiation will also 
cause the suppression of the cross section. The jet looses on average a fraction e of its 
transverse momentum when it passes through the medium, 

X (1 - e)« (2.26) 


and the jet cross section in heavy ion collisions with the creation of the medium can be 
related to the jet cross section in heavy ion collisions without the creation of the medium 
as follows. 


1 


1 

^^CNM 

1 ^ 1 


(A^bin) 

drjdpT 

(iVbin) 

Pt 

drjdpT 

prp 1 (-^bin) 

drjdpT 


(2.27) 


Here the index k = q, g labels whether the jet is quark-initiated or gluon-initiated. More 
precise calculations of the jet cross section in heavy ion collisions with the creation of the 
medium through possible factorized expressions will be investigated in the future. 

The fraction of the jet energy loss can be calculated using the medium-induced split¬ 
ting functions in SCETg- The same quantity is also studied and resummed in the vacuum 
for jets with small radii [71] 0. At leading-order, the collinear parton with momentum p 
splits into partons with momenta k and p — k. Using the anti-kx algorithm |128] with 
radius R for jet reconstruction, at this order if the angle between the final-state partons is 
smaller than R, both the particles are reconstructed in the jet and there is no lost energy. 
However, if the angle between the partons is larger than R, the more energetic parton will 
be reconstructed as the jet and the softer parton will be lost. 

^For the jet radii used in the current experiments at the LHC (typically, 0.2 < R < 0.5 which is not 
very small, and a^logi? is not large), the resummation of contributions from multiple vacuum parton 
splittings does not have a large effect in the jet energy loss calculation (Fig.4 in [71]). For medium- 
induced radiation undergoing further vacuum-like emissions, we expect a similar numerical outcome in the 
calculation of the medium-induced jet energy loss. Here we only include the lowest-order contribution, 
leaving the resummation of higher-order contributions for future work. 
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More specifically, the medium-induced energy loss of a quark jet is 


2 


rl i‘ujx{l—x)ta.n -i 

/ dxk^ + / dx{p^ — k^) / d/cj_ 

Jo Ju)x{l—x)ta.n ^ 


q^med 
' q^Q9 


{x,k_i 


+ r^Ux,k±) 


(2.28) 


whereas for gluon jets, 


OJ 


dxk^+ / dx{p^-k^) 


/ „ ^^^9 

/cja:(l—x) tan^ 


*j-)+E KiU^- *-'i) 

(2.29) 


<?>q 


Note that in the splitting g ^ gg the final-state partons are identical particles. Here R 
is the angular parameter used in the jet reconstruction, and Rq is of 0(1) in QCD which 
sets the region of the use of collinear parton splitting functions. The phase space regions 
in the above expressions are illustrated in the right panel of Fig|TJ Note the possibility of 
jet-flavor changing q ^ g oi g ^ q due to parton splitting and the restricted phase space 
constrained by the jet algorithm. The flavor-changing probability is typically small for jets 
with large transverse momenta, so is the fraction of jet energy loss. 


3. Results 

We compare our jet shape and cross section calculations to the ALICE [TI], ATLAS [121 
fT3] and CMS [TTlfTb] measurements in lead-lead collisions with nucleon-nucleon center 
of mass energy at y'sNN = 2.76 TeV. We will also make predictions for the anticipated 
measurements of both inclusive and photon-tagged jets at the y'SNN ~ 5.1 TeV LHC Run 
H. Even though the production cross section of jets associated with prompt photons is 
several orders of magnitude smaller than the inclusive jet cross section, photon-bjet events 
give a clean probe of the medium because the prompt photon can be used as a robust 
reference to study the quenching of the recoiling jet m- Also, we will see that the shape 
of photon-tagged jets will serve as a new observable to amplify the differences between the 
vacuum and the medium parton showers due to the narrow energy profile of predominantly 
quark jets that recoil against the prompt photon. Measurements of jets recoiling against a 
weak boson {W jZ') will also be interesting to perform. 

Jets were reconstructed using the anti-kx algorithm |128j with R = 0.2, 0.3, 0.4 and 0.5 
in the measurements of inclusive jet cross sections. Eor the jet shape, the CMS collaboration 
used R = 0.3. Euture measurements are planned be performed with multiple jet radii, so 
in the jet shape calculation we include both R = 0.3 and 0.5 to investigate the jet radius 
dependence. 

The inclusive jet cross sections are examined in various jet transverse momentum 
and pseudo-rapidity (ryl®*) bins. To compare with the CMS jet shape measurement, 
we impose the same cuts as follows, 

> 100 GeV , 0.3 < < 2 . (3.1) 

The region < 0.3 is excluded because of the techniques used in the background 

subtraction. In this paper the calculations for photon-bjets events are inclusive, in the 
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sense that we have allowed all possible photon kinematics consistent with the cuts on 
and Future experiments may impose cuts on the photon kinematics as well as some 
isolation criteria, and they will have to be implemented in the theory calculation. 

The differential jet shapes measured by CMS are constructed from the transverse 
momenta of the charged particles with > 1 GeV, 


Ar Nj ^ 6r 

,7=1 


(3.2) 


and the jet cone is divided into six annuli between 0 < r < 0.3 with 6r = 0.05. For the 
calculation with R = 0.5, we divide the jet cone into ten annuli between 0 < r < 0.5 with 
the same 5r. 

It should be noted that, the ATLAS measurements used calorimeter towers with Arj x 
A(p = 0.1 X 0.1 for jet reconstructions, where the detector granularity may start to play 
a role for the reconstruction of ii = 0.2 jets. Also, ALICE measured charged jets and 
CMS measured the jet shape using charged particles. The above details are not taken into 
account in the theory calculation yet, however, the difference will be examined in the future 
as the theory and the experiments progress. 

For the cross section calculations of both inclusive and photon-tagged jets, we use 
the CTEQ5M parton distribution functions (PDEs) [129j and the leading-order 0(0^) and 
0{asaern) QCD partonic cross section results. The cross sections in Pb-|-Pb collisions are 
calculated by implementing the initial-state cold nuclear matter effects and the final-state 
jet energy loss caused by the medium-induced splitting. To examine the CNM effect and 
its theoretical uncertainty, we consider three different cases in |122j with the character¬ 
istic RMS medium momentum exchange scales //cnm = 0 GeV, pcNM = 0.18 GeV and 
/7cNM = 0.35 GeV corresponding to different sizes of initial-state energy loss effects. Eor 
the consistency of the theory calculations of the jet-medium interaction, we use the same 
splitting kernels which were used in the description of inclusive particle production and 
the detailed expressions can be found in [llOj . The coupling p |§ between the jet and the 
medium is typically chosen to be around g = 2.0, and we estimate the theory uncertainty 
using the window g = 2.0 ± 0.2 in most of the cases. 

As we will see below, varying the medium momentum exchange scale pcNM between 
0 GeV and 0.35 GeV has a large effect on the nuclear modification factor Raa (defined in 
Eq. (I1.3p l at high pT- Also, an 0(10%) variation of the coupling g can cause considerable 
change in the Raa values. These suggest that the inclusive jet cross section is sensitive 
to both initial state and final state effects. However, later we will see that the jet shape 
is less sensitive to the variation of pcNM therefore allowing us to precisely probe the final 
state, jet-medium interaction. This in turn can allow us to disentangle and more reliably 
constrain the CNM effect using the Raa data. The typical value of g = 2.0 used in the 
calculation is consistent with the jet shape modification measurement. It should also be 

®as = <7^/47r is the strong coupling constant associated with the jet-medium interaction. We estimate 
the uncertainty by varying the interaction strength with fixed g’s within g = 2.0±0.2. The detailed medium 
effects absorbed in the coupling g will be investigated in the future. 
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Figure 2: Comparison of theoretical calculations for the nuclear modification factor Raa of in¬ 
clusive jets as a function of the jet transverse momentum to experimental data in central Pb-|-Pb 
collisions at ^snn = 2.76 TeV at the LHC. Bands correspond to the theoretical uncertainty esti¬ 
mated by varying the coupling between the jet and the medium {g = 2.0 ± 0.2). The blue band 
corresponds to the calculations with p.cNM = 0 GeV, the green band corresponds to the calculations 
with ficNM = 0.18 GeV and the red band corresponds to the calculations with /tcnm = 0.35 GeV. 
Left panel: comparison to the ATLAS measurement [13] with R = 0.4. Right panel: comparison to 
the CMS preliminary result m with R = 0.3. 


noted that, in this paper we focus on the effect of radiative energy loss without including 
the effect of collisional energy loss in the jet-medium interaction. While collisional energy 
loss may not affect the jet shape much, it can give significant contributions to the Raa 
suppression. Detailed studies of the CNM effect and the collisional energy loss will be 
performed in the future. 

The jet shape calculated in SCET are averaged with the jet cross section to produce the 
final jet shape to be compared with the measurements. In the calculations of the differential 
jet shapes, we include the one-loop jet energy functions and their medium modifications. 
For the renormalization group evolution we include the two-loop cusp anomalous dimen¬ 
sions (To and Ti) and the one-loop anomalous dimensions of the jet energy functions, 
as well as the two-loop running of the strong coupling constant with as{mz) = 0.1172 |64j . 
This allows us to resum the jet shape at next-to-leading logarithmic accuracy (NLL). The 
theoretical uncertainty will be estimated conventionally by varying the jet energy scales 
within a factor of 2. 

3.1 Jet shapes and cross sections in Pb-|-Pb collisions at .^snn = 2.76 TeV 

We first present theoretical calculations and the comparisons to the measurements in 
Pb-|-Pb collisions at .^snn = 2.76 TeV. Fig. [2] shows the results of the nuclear modifi¬ 
cation factor Raa of inclusive jet productions for jets with pseudo-rapidity |r/| < 2 in the 
0 - 10% centrality class as a function of jet transverse momenta pT- We choose to present 
the calculation for jets with pT > 50 GeV. For jets with small pT, hadronization effects 
can play a significant role and our perturbative calculations will have to include these 
non-perturbative contributions. On the other hand, with the huge underlying event back¬ 
grounds, the reconstruction of lower pT jets in experiments becomes increasingly difficult. 
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Figure 3: Comparison of theoretical calculations for the suppression of inclusive jets with R = 0.3 
as a function of the jet transverse momentum to experimental data in Pb+Pb collisions at .^/snn = 
2.76 TeV. The curves are the jet Raa calculations for central colliisons, but we have included 
the data from the ALICE charged jets Rcp [H], ATLAS calorimeter jet Rqp [H] and CMS jet 
Raa na measurements for comparison, assuming that the suppression in peripheral collisions is 
small. Shown are the same studies as in Fig. 2 of the theory sensitivity on the coupling g between 
the jet and the QGP medium, as well as the momentum transfer ^cnm between the initial-state 
parton and the nucleus. 

Fig. [2] and E] present the theoretical sensitivity of Raa on cold nuclear matter effects 
caused by initial state parton energy loss. The magnitude of CNM effects depends strongly 
on the magnitude of the momentum transfer between the incident partons and the nucleus. 
We show three different cases with /icNM = 0 GeV (represented by the blue bands), pcNM = 
0.18 GeV (represented by the green bands), and /icNM = 0.35 GeV (represented by the red 
bands). Here, pcNM = which is the root-mean-square of the transverse momentum 

transfer per scattering between the incident hard parton and the constituents of the nucleus. 
On the other hand, the strength of the final state jet-medium interaction is controlled by 
the coupling g between the collinear partons and the QGP medium. The width of the 
bands in the plots corresponds to the variation of g = 2.0 ± 0.2. 

In the left panel of Fig. [2l we present the comparison of the theoretical results to 
the ATLAS Raa measurements with R = 0.4. We find that the final state jet-medium 
interaction alone is not sufficient to describe the suppression of the cross section at high 
transverse momenta, and the inclusion of CNM effects with /icNM = 0.35 GeV provides a 
good description of the experimental data (red band). The comparison of our calculations 
to the preliminary CMS data with R = 0.3 is shown in the right panel of Fig. [5] with similar 
conclusions. 

In Fig. [3l we include a comparison to the ALICE Rcp measurement of the suppression 
of R = 0.3 charged jets HH, i.e. jets reconstructed with the charged particle tracks. We 
find that the calculation with pcNM = 0.35 GeV is consistent with the quenching of charged 
jets in the comparatively low pp region. In the high intra-jet particle multiplicity limit, 
from isospin symmetry we expect that charged jets will exhibit similar characteristics as 
the full reconstructed calorimeter jets, therefore with similar cross section suppression. A 
comparison to the ATLAS Rcp measurement with R = 0.3 is also presented in the same 
plot. 
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Figure 4: Left panel: comparison of theoretical calculations for the nuclear modification factor 
Raa of inclusive jets as a function of the jet transverse momentum to experimental data in y^SNN = 
2.76 TeV Pb+Pb collisions at the LHC, with different collision centralities. Two centrality classes, 0 
- 10% and 30 - 40%, are considered. The bands corresponds to the variation of the coupling g = 2± 
0.2 between the jet and the medium in the calculations, and CNM effects with /xcnm = 0.35 GeV are 
implemented. The data is from ATLAS m with R = 0.4. Right panel: comparison of theoretical 
calculations for the ratios of jet Raa’s within different pseudo-rapidity bins as a function of the 
jet transverse momentum to the ATLAS experimental data m in central Pb-|-Pb collisions. The 
blue band corresponds to the ratio Raa{ 0-0 < Ivl < 0.8)/i?AA(0.0 < \r]\ < 2.0), and the red band 
corresponds to the ratio Raa{ 0-& < \v\ < 2.0)/i?AA(0.0 < \r]\ < 2.0). The derived data corresponds 
to i?AA(0.3 < |? 7 | < 0.8)/i?AA(0.0 < \r]\ < 2.1) and i?AA(l-2 < \ri\ < 2.1)/i?^^(0.0 < \r]\ < 2.1). 


Next, we examine the centrality and the rapidity dependence of inclusive jet suppres¬ 
sions in Pb-|-Pb collisions at y^snn = 2.76 TeV. In the left panel of Fig. 01 we show the 
theoretical calculations of the inclusive jet Raa in central and mid-peripheral collisions and 
the comparison to the ATLAS measurements in the 0 - 10% and 30 - 40% centrality classes. 
The framework for the evaluation of the medium-induced splitting functions and their im¬ 
plementation in the calculation of the jet cross sections in heavy ion collisions capture 
the relative contributions of the vacuum and the in-medium parton showers as a function 
of the size and the density/temperature of the QGP, which is controlled by the collision 
centrality. The attenuation of jet production in peripheral collisions is smaller than the 
attenuation in central collisions, which is consistent with the one seen in the quenching of 
inclusive hadrons m- The right panel illustrates the dependence of jet suppression on 
the pseudo-rapidity r] of jets. Shown is the Raa of inclusive jets as a function of px in 
the 0 < \ri\ < 0.8 and 0.8 < |t/| < 2.0 intervals relative to the Raa in the 0 < |7/| < 2.1 
interval. Note that within experimental uncertainties these ratios are consistent with unity. 
Few percent differences in the theoretical calculations capture the qualitative feature in the 
ATLAS data. 

An important check of the calculations is to examine the jet radius R dependence of 
the jet cross section and its medium suppression. We illustrate the different suppression 
patterns with different jet radii in the left panel of Fig. [5l We consider central Pb-|-Pb 
collisions at the LHC for illustration. The smaller range of the coupling g = 2.0 ±0.1 is 
used to estimate the theory uncertainty for better separation of the results corresponding 
to different radii. The red band represents R = 0.2, the orange band R = 0.3, the green 
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Figure 5: Left panel: theoretical calculations for the nuclear modification factor Raa of inclusive 
jets with R = 0.2, 0.3, 0.4 and 0.5 as a function of the jet transverse momentum in central Pb+Pb 
collisions at .^snn = 2.76 TeV at the LHC. The calculations implement CNM effects with p.cNM = 
0.35 GeV, and a reduced range of the coupling between the jet and the medium (5 = 2 ± 0.1) is 
used for estimating theoretical uncertainties. This way the bands do not overlap much and the jet- 
radius dependence of jet quenching can be better illustrated. Right panel: comparison of theoretical 
calculations for the central-to-peripheral Rcp ratios as a function of the jet transverse momentum 
for inclusive jets with different jet radii to data in central Pb-fPb collisions at -s/s^ = 2.76 TeV 
at the LHC. The bands and the data correspond to Rcp{R = 0.3, 0.4, 0.5)/Rcp{R = 0.2). The 
data is from ATLAS [T^ . 


band R = 0.4 and the blue band R = 0.5. The smaller the jet radius, the larger the 
suppression of the cross section. The jet radius dependence of Raa is indeed observed by 
the ATLAS collaboration [12] and shown in the right panel of Fig. [5j More specifically, the 
central-to-peripheral ratio Rcp of jet cross sections as a function of pp, dehned as 


R^pipp) = 


d(^AAiPT,R) 

dr](Ppp 



d(T'AAiPT,R) 

dpcPpp 


RTa{pt,R) 
RTa(pt,R) ■ 


(3.3) 


is compared to the one for R = 0.2. In Eq. (|3.3I1 (A'bf”') and are the mean of the 

number of binary nucleon-nucleon scattering in central and peripheral Pb-|-Pb collisions. 
The calculation predicts qualitatively the transverse momentum dependence of the Rcp 
ratios and provides a quantitative description of this observable when the two radii are 
sufficiently different. For small radii the calculation over-predicts the difference in the 
quenching patterns of inclusive jets. Such discrepancy may be caused by the fact that we 
only evaluate the medium-induced splitting functions at the lowest non-trivial order, or 
that we neglect further dissipation of the jet energy through collisional processes in the 
medium The resummation of log R for small-i? jet cross section can also play a 

role mills]. New experimental measurements will be very useful to further examine this 
jet radius dependence. 

The jet shape can give complimentary information of the in-medium parton shower 
beyond the study of the jet cross section. Since the differential jet shape is normalized 
by p{r)dr = 1, we expect any enhancement (attenuation) at r Ri i? to be correlated 
with the attenuation (enhancement) at small/intermediate values of r. This behavior 
was qualitatively seen in the attempt to calculate the jet shape modification using the 
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Figure 6: Comparison of theoretical calculations for the modification of differential jet shapes 
of inclusive jets in Pb+Pb central (left panel) and peripheral (right panel) collisions at .^/snn = 
2.76 TeV at the LHC. The modification is presented as the ratio of the jet shapes (r)/{r) 
in Pb-Pb and p-p collisions. We impose cuts on the jet transverse momentum px > 100 GeV and 
pseudo-rapidity 0.3 < \r]\ < 2.0 of jets. The coupling between the jet and the medium is fixed at 
g = 2. The blue band corresponds to the calculations including only the CNM effects, the red band 
adds the jet-medium interaction but with the jet-by-jet shape modification turned off, and the green 
band correspond to the full calculation. The theoretical uncertainty is estimated by varying the jet 
energy scales < p < 2pj^ in the calculations. The data is from CMS [2] with R = 0.3. 


traditional energy loss approach [39], but the theory was not able to provide quantitative 
description of the CMS data HD. 

In Fig. E] we compare the theoretical calculations of the jet shape modification in 
central and mid-peripheral Pb-|-Pb collisions at = 2.76 TeV at the LHC to the CMS 

data m with R = 0.3. In the calculations shown here, we fix the coupling between the 
collinear partons and the medium to be <7 = 2 and estimate the theory uncertainty from 
the variation of the jet energy scales. Jet shape modifications can partly come from the 
change of the relative fraction of quark jets and gluon jets since quark and gluon jet shapes 
differ. The blue band represents the result with only the CNM effect, and the deviation 
from 1 is very small. This implies that the CNM effect does not significantly affect the 
relative fractions of quark and gluon jets therefore the jet shape does not change much, 
whereas it can result in large Raa suppressions in the absolute inclusive jet cross sections. 
We show the blue band only for central collisions to demonstrate this point, and the jet 
shape is not as sensitive to the CNM effect implemented in this paper as the inclusive jet 
cross section is. The red band corresponds to the calculation including the cross section 
suppression but assuming that the jet-by-jet shape remains the same as in proton collisions. 
Here, the fraction of quark jets significantly increases since the cross section of gluon jets 
are more suppressed by the jet-medium interaction. This leads to the narrowing of the jet 
energy profile. The other contribution to jet shape modifications comes directly from the 
modification of the in-medium parton shower which leads to the broadening of jets. The 
green band includes all the above physics inputs which result in the attenuation at mid r 
and the enhancement at r ~ i? of the jet shape. This gives, for the first time, a quantitative 
understanding of the jet shape modification in heavy ion reactions. The calculations for 
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Figure 7: Theoretical predictions of the nuclear modification factor Raa as a function of jet 
transverse momentum pT for inclusive jets, with R — 0.2 (left panel) and R = 0.4 (right panel) in 
Pb+Pb collisions at ^snn ~ 5.1 TeV at the LHC. The coupling between the jet and the medium 
is varied to estimate the theoretical uncertainties {g = 2.0 ± 0.2), and CNM effects with /xcnm = 
0.35 GeV are implemented in the calculations. The red band corresponds to the collision centrality 
0 - 10% and the blue band corresponds to the collision centrality 30 - 50%. 

central collisions are shown in the left panel of Fig. [6l and the ones for peripheral collisions 
are shown in the right panel. 

3.2 Theoretical predictions for jet shapes and cross sections in Pb+Pb collisions 

at +snn ~ 5.1 TeV 

In this section we make predictions for the anticipated jet shape and cross section mea¬ 
surements at the +snn ~ 5.1 TeV LHC Run II. We show the predictions for both inclusive 
and photon-tagged jets and compare the two. As we will see, by studying processes with 
different compositions of quark jets and gluon jets, we are able to examine the physical 
understanding of jet shapes and cross sections which allows us to use quark and gluon jets 
as independent probes of the QGP. 

Fig. [7]shows the predictions of the nuclear modification factor Raa of inclusive jet cross 
sections in central (red band, to be compared with future measurements with centrality 
0 - 10%) and mid-peripheral (blue band, to be compared with future measurements with 
centrality 30 - 50%) collisions at +snn ~ 5.1 TeV. The left panel corresponds to jets with 
R = 0.2 and the right panel shows the R = 0.4 results. Larger suppression is seen in 
central collisions and for jets with smaller radius. Again the theoretical uncertainties are 
estimated by varying the coupling g = 2.0 ± 0.2. All the plots in this section include the 
implementation of CNM effects with gcNM = 0.35 GeV since the cases of gcNM = 0 GeV 
and gcNM = 0.18 GeV are not consistent with the +snn = 2.76 TeV data. 

Fig. [8]shows the predictions of the nuclear modification factor Raa of jet cross sections 
for photon+jet events in central (orange band) and mid-peripheral (purple band) collisions 
at +snn ~ 5.1 TeV. The theoretical uncertainty is estimated with g = 2.0 ± 0.2. The jet 
suppression shows similar dependence on the centrality and the jet radius as the one for 
inclusive jets. Note that, the cross section suppression as a function of pp is more flat for 
photon-tagged jets compared to the one for inclusive jets. This can be seen more clearly 
in Fig. Oin which we superimpose Fig. [7] and Fig. [Altogether. Also, for each choice of jet 
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Figure 8: Theoretical predictions of the nuclear modification factor Raa as a function of jet 
transverse momentum pT for jets recoiling against a prompt photon, with R — 0.2 (left panel) 
and R = 0.4 (right panel) in Pb+Pb collisions at .^snn ~ 5.1 TeV at the LHC. The theoretical 
uncertainty estimation and the implementation of CNM effects are as in Figure 7. The orange band 
corresponds to the collision centrality 0 - 10% and the purple band corresponds to the collision 
centrality 30 - 50%. 
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Figure 9: Comparison of the theoretical predictions of the nuclear modification factor Raa as a 
function of jet transverse momentum px for inclusive jets and photon-tagged jets, with R = 0.2 (left 
panel) and R = 0.4 (right panel) in central and mid-peripheral Pb-|-Pb collisions at y^SNN ~ 5.1 
TeV at the LHC. The coupling between the jet and the QCD medium is between g = 2.0 ±0.2, and 
CNM effects with /xcnm = 0.35 GeV are implemented. The bands correspond to: inclusive jets, 0 - 
10% (red), inclusive jets, 30 - 50% (blue), photon ± jet, 0 - 10% (orange), photon ± jet, 30 - 50% 
(purple). 


radius and collision centrality, the Raa^s for inclusive jets and photon-tagged jets converge 
at high pp. This comes from the fact that, jets with higher pp will be less modified by 
the interaction with the medium. In the high pp limit the effect from final-state jet energy 
loss will decrease to zero because it is expected to be suppressed by the LPM effect and 
parametrically proportional to the ratio between the medium temperature and the jet pp. 
In this regime the cross section suppression is mainly due to the initial-state CNM effects 
which we use the same implementations in this section. At lower pp, the cross section 
suppression for photon-tagged jets is smaller than the suppression for inclusive jets. This 
is because the jets in photon+jet events are predominantly quark jets and the medium 
causes less jet energy loss as opposed to gluon jets. 
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Figure 10: Comparison of the theoretical predictions of the differential jet shapes for inclusive 
jets (blue bands) and photon-tagged jets (red bands), with R — 0.3 (left panel) and R —0.5 (right 
panel) in y/s Ri 5.1 TeV p-l-p collisions, with next-to-leading logarithmic accuracy. We impose cuts 
on the jet transverse momentum px > 100 GeV and pseudo-rapidity 0.3 < \ri\ < 2. The bands 
correspond to the theoretical uncertainties estimated by varying the jet energy scales between 
< fi < in the calculations. Note the narrower energy profile of photon-tagged jets which 
are mostly quark-initiated. 


We can see more clearly the difference of quark-jet and gluon-jet fractions for inclusive 
jets and photon-tagged jets in Fig. [101 The figure shows the jet shapes for inclusive jets 
(blue band) and photon-tagged jets (red band) in proton collisions at ^/s « 5.1 TeV at the 
LHC. Here we impose the same cuts on the jet transverse momentum pp > 100 GeV and 
pseudo-rapidity 0.3 < \r]\ < 2.0 as the ones in the CMS measurements at .^snn = 2.76 TeV. 
The left panel shows the jet shape for jets with R = 0.3 and the right panel shows the one 
for R = 0.5. In either case, we see that the shapes of photon-tagged jets are clearly narrower 
than the shapes of inclusive jets, which indicates that photon-tagged jets are predominately 
quark-initiated, and inclusive jets have a significant fraction of gluon-initiated ones. 

With the predictions for the baseline jet shapes in proton collisions at ^/s ~ 5.1 TeV, 
hnally we now present the predictions for the modification of jet shapes in Pb-|-Pb collisions 
at y/sNN ~ 5.1 TeV. Fig. [TT] shows the jet shape modification for both inclusive jets 
(green band) and photon-tagged jets (blue band) with R = 0.3 in central (left panel) and 
mid-peripheral (right panel) collisions. The theoretical uncertainty is again estimated by 
varying the jet energy scales in the SCET calculations. While the modification of inclusive 
jet shapes at .^/snn ~ 5.1 TeV is similar to the modification at y^snn = 2.76 TeV, we 
see that the modihcation of photon-tagged jet shapes has a different pattern, which shows 
clearly the broadening of jets. This is consistent with the fact that photon-tagged jets are 
predominately quark-initiated so that different suppressions of quark and gluon jet cross 
sections in the medium do not modify the jet shape as much: the gluon-jet contribution to 
the photon-tagged jet shape is not much. 

Fig. [12] shows similarly the predictions for the jet shape modification of inclusive jets 
and photon-tagged jets with i? = 0.5 at the « 5.1 TeV LHC. The features discussed 

in Fig. [TT] remain the same. However, because of the larger jet radius the jet cross section 
suppression is smaller for R = 0.5 jets. This leads to the clearer effect of jet broadening. 
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Figure 11: Theoretical predictions of the modification of differential jet shapes for inclusive jets 
(green band) and photon-tagged jets (blue band), with R = 0.3 in central (left panel) and mid¬ 
peripheral (right panel) Pb-|-Pb collisions at ^snn ~ 5.1 TeV at the LHC. We impose cuts on 
the jet transverse momentum pT > 100 GeV and pseudo-rapidity 0.3 < |r;| < 2.0. The coupling 
between the jet and the medium is fixed at g = 2 . 0 , and the bands correspond to the theoretical 
uncertainties estimated by varying the jet energy scales. 
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Figure 12: Theoretical predictions of the modification of differential jet shapes for inclusive jets 
(green band) and photon-tagged jets (blue band), with R = 0.5 in central (left panel) and mid¬ 
peripheral (right panel) Pb+Pb collisions at ysNN ~ 5.1 TeV at the LHC. The cuts ofpr > 100 GeV 
and 0.3 < I 77 I < 2.0 are imposed, and the bands represent the theoretical uncertainties in the 
calculations. 


4. Summary and discussion 

In this paper, using soft collinear effective theory and its extension SCETq to include 
collinear parton interactions in dense QCD matter via Glauber gluon exchange, we calcu¬ 
late the modification of jet cross sections and differential jet shapes in lead-lead collisions 
at the LHC. The medium contributions to the jet energy function are evaluated using the 
medium-induced splitting kernels obtained in SCETq. This work emphasizes the consis¬ 
tent theoretical descriptions of hadron and jet observables in heavy ion collisions. The 
results presented in this paper also include cold nuclear matter (CNM) effects, which have 
been shown to affect the nuclear modification factor Raa and seen more clearly at high 
transverse momentum. We find that the calculations can describe well the centrality and 
the Pt dependence of the inclusive jet suppression observed by the ALICE, ATLAS and 
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CMS collaborations in y^snn = 2.76 TeV Pb+Pb collisions at the LHC, with the cou¬ 
pling between the jet and the medium g = 2.0 ± 0.2 and CNM effects with gcNM = 0.35 
GeV. The jet radius R dependence of jet quenching is also qualitatively consistent with 
the measurements. 

We also calculate the jet shape in nucleus-nucleus collisions. We find that the non¬ 
trivial behavior of the jet shape modification is caused by both the different quark and 
gluon jet cross section suppressions and the jet-by-jet broadening. The cross section of 
gluon-initiated jets is more suppressed, which enhances the fraction of quark-initiated jets 
having a narrower energy profile. This causes the attenuation of the jet shape in the 
mid r region. On the other hand, the broadening of jets results in the enhancement of 
the jet shape near the periphery of the jet. The calculation provides for the first time a 
quantitative description of the jet shape modification in Pb-|-Pb collisions at the LHC. We 
would also like to emphasize that, the jet shape and substructure observables in general 
are less sensitive to initial state effects. This allows us to cleanly disentangle final state 
jet-medium interactions and directly probes the medium properties. A cleaner study of 
the initial state effects in A-|-A collisions can then be done. 

For the upcoming LHC Run H measurements, we present theoretical predictions of 
the jet shape modification and the cross section suppression for inclusive jets and photon- 
tagged jets. We find that the cross section suppression at high pT can provide information 
about the cold nuclear matter effects. Since photon-tagged jets are predominately quark- 
initiated, the cross section is expected to be less suppressed compared to inclusive jets. On 
the other hand, the broadening of the photon-tagged jet is more apparent. 

With the understanding of the transverse distribution of intra-jet particles through the 
studies of jet shapes and cross sections, it would be interesting to study the modification of 
their longitudinal distribution to gain orthogonal information about the in-medium parton 
shower |17[ll8[fT3Tlll32| . We leave the studies of jet fragmentation function modification 
using the SCET formulation |12n[ll33[ll34j . as well as the modification of jet masses |135L 
I138j which probes the jet formation mechanism at the soft scale, for future work. 
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